NONLOCAL CONDENSATES AND CURRENT-CURRENT CORRELATORS 
WITHIN THE INSTANTON LIQUID MODEL 



A.E. DOROKHOV 

The quark and gluon nonlocal condensates and current-current correlators are discussed within the 
instanton liquid model. 

O 

O ' 1 Gluon field strength correlator. 

(N ' 

■ The non-perturbative vacuum of QCD is densely populated by long - wave fluctuations 

of gluon and quark fields. The order parameters of this complicated state are character- 
ized by the vacuum matrix elements of various singlet combinations of quark and gluon 

fields, condensates: (: qq :), ^: F^ U F^ U :W: q{crp V F^ u ^)q :^ ,etc. The nonzero quark 

condensate (: qq :) is responsible for the spontaneous breakdown of chiral symmetry, 
and its value was estimated a long time ago within the current algebra approach. The 

nonzero gluon condensate ^: F^ U F^ U through trace anomaly provides the mass scale 

for hadrons. Its value was estimated within the QCD sum rule (SR) approach^ and first 

evidence on its existence has been obtained irJ21. The values of low - dimensional con- 
densates were obtained phenomenologically from the QCD SR analysis of the current- 
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current correlators in various hadron channels. In it was proposed to study the effects 
of the gluon condensate on a quarkonium state in order to estimate its absolute value. 

fH If |/r| 

The nonlocal vacuum condensates or vacuum correlators describe the distribution 



of quarks and gluons in the non-perturbative vacuum. Physically, it means that vacuum 
quarks and gluons can flow through the vacuum with nonzero momentum. From this 
point of view the standard vacuum expectation values (VEVs) like (: qq :), (: qD 2 q :) , 

(: g 2 F 2 :) , ... appear as expansion coefficients of the quark M(x) = I : q(0)E(0, x)q(x) 



and gluon D^ v,pa {x) correlators in a Taylor series in the variable x 2 /4. The correlator 
j)nv,pcr ^ Q f gi uon i c field strengths 

DW{x -y) = (: TrF^(x)E(x, y)F^(y)E(y, x) :) , (1) 

may be parameterized in the form consistent with general requirements of the gauge and 
Lorentz symmetries a 

D^{x) = -L(-.F 2 :){(5p P 5 UCT -5 IMCT 5 up )[D(x 2 ) + D 1 (x 2 )]+ (2) 
, dD 1 (x 2 ) 

+ \X iiX pO va X^XfjOyp + X v X a Opp X V X pO p a ) ~Qx^ 

where E(x,y) = Pexp(i f£ A fJi (z)dz fJ ') is the path-ordered Schwinger phase factor 
(the integration is performed along the straight line) required for gauge invariance and 

Ap{z) = A^z)^, Fp U (x) = F^(x)^-, F- U (x) = dpA-(x)-d u A-(x)+f abc A^(x)AUx). 



The P— exponential ensures the parallel transport of color from one point to other. In 

©, (: F 2 :) = F^(0)i^(0) :\ is a gluon condensate, and D(x 2 ) and D^x 2 ) are 

invariant functions which characterize nonlocal properties of the condensate in differ- 
ent projections. The form factors are normalized at zero by the conditions D(0) = k, 
Di(0) = 1 — k, that depend on the dynamics considered. For example, for the self- 
dual fields k = 1, while in the Abelian theory without monopoles the Bianchi identity 
provides k = 0. 

one has shown that the instanton model of the QCD vacuum provides a way 
to construct nonlocal vacuum condensates. Within the effective single instanton (SI) 
approximation one has obtained the expressions for the nonlocal gluon Dj U,peT (x) and 
quark Mr(x) condensates and derived the average virtualities of quarks A 2 and gluons 
X 2 in the QCD vacuum. The behavior of the correlation functions demonstrates that in 
the SI approximation the model of nonlocal condensates can well reproduce the behavior 
of the quark and gluon correlators at short distances. 

InM^l, it was suggested that the instanton A^ 1 (x) is developed in the physical vacuum 
field bn(x) interpolating large-scale vacuum fluctuations. One has found that at small 
distances the instanton field dominates, and at large distances it decreases exponentially 

- (r/g |x|) 3 / 2 unlike the power-like decreasing SI. It is important to note that 
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the form of this asymptotics is independent of the model for the background field and the 
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driven parameter r\ g ~ y - ^ 2 -^R (F b j^j , where R is the correlation length and 

{F 2 ) b is the background field contribution to the gluon condensate, only weakly depends 

on it. This solution is called the constrained instanton (CI)^^ 

The knowledge of the constraint-independent parts of CI allowed us to construct the 
solution in the ansatz form 

= Vt^ 9 (* 2 ) , V 9 (* 2 ) = (3) 

where the notation p 2 (x 2 ) = a i / 3 rj 2 x 2 K 4: / 3 ^(r) g x) 3 / 2 , p 2 (0) = p 2 is introduced. By 

translational invariance the center of CI can be shifted in © from the origin to an arbitrary 
position xq\ x — > x — xq. The constrained instanton model introduces two characteristic 
scales (correlation lengths). The short distance behavior of the correlation functions is 
proportional to the instanton size and dominated by the single instanton contribution. The 
large correlation length is physically related to the confinement size R. One sees that the 
D (x 2 ) structure is close to the SI induced function with the exponential asymptotics 
being developed at large distances. 

The gluon field strength correlation functions D (x 2 ) and D\ {x 2 ) have been found 
i n QSl an( j t ^ e resu it s are shown in Fig. ^ We see that the results of the constrained 
instanton model calculations are in a remarkable agreement with the results of lattice 
QCD simulations: 1) The D% (x 2 ) structure is much smaller than the D (x 2 ) function 

| 12 | 10 [ 2) jj, e correlation length of the gluon field correlator is much smaller than for the 
quark field correlator! ^ | 8[ 3) t he pa th depe ndence of the gluon field strength correlator 
has similar features in both approach 14 | 15 | 



Within the constrained instanton model there is a firm predictions^, that if the local- 
ized gluon classical field (constrained instanton) of the general form 

A^ a (x) = r(t^ a {x\ (4) 

with properties <p g (x 2 — ► 0) = 1 and tp g (x 2 — > oo) — > faster than any power of 1/x, 
then the correlation functions must have zeros in preasymptotic region on the background 
of their exponential decay. Thus the predicted asymptotics is of kind 

D(x) ~ exp (— Ax) cos (bx + S(x)) (5) 

rather than pure exponential one. The appearance of oscillations is possible exclusively 
due to the Schwinger string phase factor (link) and so is very interesting 51 . For inter- 
pretation of these oscillations as specific "confining" behaviour of the propagator-like 

functions seeQS. 
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Figure 1: The amplitudes D (top lines) and Di (bottom 
lines) (all normalized by -D(O)) versus physical distance 
x, for the instanton size p = 0.3 fm and parameters 
(PVg) 2 = (solid lines) and (prj g ) 2 = 1 (dashed lines). 
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Figure 2: The Adler function from the ILM contri- 
butions: dynamical quark loop (short dashed), quark 
+ chiral loops + vector mesons (full line) versus the 
ALEPH data (dashed). The dash-dotted line is the pre- 
diction of the constituent quark model (extended NJL) 
and the dotted line is the asymptotic freedom prediction, 

1/47T 2 . 
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The instanton model predicts the behaviour of nonperturbative part of gluon correla- 
tion functions in the short and intermediate region assuming that it is dominated by in- 
stanton vacuum component, while the large-scale asymptotics is dominated by the back- 
ground field. 



a In many works devoted to theoretical analysis of the gluon field strength correlators the Schwinger link 
is neglected and so the important effect of oscillations is missed. 



The gluon correlators are the base elements of the stochastic model of vacuum El. 
They are also oft en used in the description of high-energy hadron diffractive scattering 

amplitude J^ I ^ I. However, as it was shown irU^the reduction of these amplitudes to the 

gluon correlators is only valid for the Abelian model like that considered irU-^. As explicit 
calculations made within the instanton model show in general case these amplitudes are 
process dependentPSESI 



2 Hadronic current-current correlators. 

The transition from perturbative regime of QCD to nonperturbative one has yet remained 
under discussion. At high momenta the fundamental degrees of freedom are almost mass- 
less quarks. At low momenta the nonperturbative regime is adequately described in terms 
of constituent quarks with masses dynamically generated by spontaneous breaking of chi- 
ral symmetry. The instanton model of QCD vacuum provides the mechanism of dynam- 
ical quark dressing in the background of instanton vacuum and leads to generation of the 
momentum dependent quark mass that interpolates these two extremes. Still it is not clear 
how an intuitive picture of this transition may be tested at the level of observables. Be- 
low we demonstrate that the Adler function and the amplitudes related to the anomalous 
triangle diagram depending on spacelike momenta may serve as the appropriate quantity. 

The Adler function defined as the logarithmic derivative of the current-current corre- 
lator can be extracted from the experimental data of ALEPH and OPAL collaborations 
on inclusive hadronic r decays. From theoretical point of view it is well known that in 
high-energy asymptotically free limit the Adler function calculated for massless quarks 
is a nonzero constant. From the other side in the constituent quark model this function is 
zero at zero virtuality. Thus the transition of the Adler function from its constant asymp- 
totic behaviour to zero is very indicative concerning the nontrivial QCD dynamics at in- 
termediate momenta. Below we intend to show that the instanton liquid model which is 

RTI 

a nonlocal chiral quark model (ILM) describes this transition correctly . The use in the 
calculations of a covariant nonlocal low-energy quark model based on the self-consistent 
approach to the dynamics of quarks has many attractive features as it preserves the gauge 
invariance, is consistent with the low-energy theorems, as well as takes into account the 
large-distance dynamics controlled by the hadronic bound states. 

In ILM in the chiral limit the vector currents correlator has a transverse character^! 

K» (Q 2 ) = - ?f) nr (q 2 ) , (6) 

where the polarization function is given by the sum of the dynamical quark loop, the 
intermediate vector mesons and the higher order mesonic loop contributions. The lowest 
order spectral representation of the polarization function consists of zero width vector 
resonances and two-meson states. The dynamical quark loop under condition of analyti- 
cal confinement has no singularities in physical space of momenta. 

The dominant contribution to the vector current correlator at space-like momentum 
transfer is given by the loop of light quarks with dynamical momentum dependent mass 
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8iV c f d 4 k M(k) 



Q 2 J (2vr) 4 D (k) 



M'(k) - -k 2 ± M {2) (k,k + Q,k) 



where the notations k± = k±Q/2, k\ = k+k-- {k+q) ^- q \ M± = M(k±), D{k) = 

k 2 + M 2 (k), D± = D(k±) are used. We also introduce the finite-difference derivatives 
defined for an arbitrary function F (k) as 



FW(k,k') 



F(k') - F(k) 
k' 2 - k 2 ■ 



The expression for ny Loop (Q 2 ) is formally divergent and needs proper regularization 
and renormalization procedures which are symbolically noted by [■■] ren for the divergent 
term. At the same time the corresponding Adler function is well defined and finite. 

Also we have checked that there is no pole in the vector correlator as Q 2 — ► 0, which 
simply means that photon remains massless with inclusion of strong interaction. In the 
limiting cases the corresponding Adler function satisfies general requirements of QCD 
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The leading high Q 2 asymptotics comes from the 
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term in Q, while 



the subleading asymptotics is driven by "tachionic" or < A 2 > term with coefficientPZ 



or 
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(10) 



It is possible to integrate Eq. dTOb in the dilute liquid approximation, u » M 2 (u), 

N r 
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GeV z , 



(11) 



which is close to exact result|22 and phenomenological estimate fromE3. 
By using set of parameters found in 

ILMEHthe 

Adler function in the vector channel 

is presented in Fig. 121 



Within the context of ILM, the integrals over the momentum are calculated by transforming the integra- 



tion variables into the Euclidean space, (fc° -> ik 4 , k 2 



-k 2 ). 



3 VAV correlator 



Since discovery of anomalous properties^EIlof the triangle diagram with incoming two 



vector and one axial- vector currents many new interesting results have been gained. 
Recently the interest to triangle diagram has been renewed due to the problem of accurate 
calculation of higher order hadronic contributions to muon anomalous magnetic moment 
via the light-by-light scattering process, that cannot be expressed as a convolution of 
experimentally accessible observables and need to be estimated from theory. 

The triangle amplitude involving the axial current A and two electromagnetic currents 
(one soft V and one virtual V), which can be viewed as a mixing between the axial 
and vecto r currents in the external electromagnetic field, were considered recently in 
| 27 | 28 | 29 | jhis amplitude can be written as a correlator of the axial current jjj and two 
vector currents j v and jy. 

f^x = ~J d 4 *dV^-^ <0| T{j v (x)Uy)j 5 x (0)}\0) , (12) 

with the tilted current being for the soft momentum photon vertex. In the specific kine- 
matics when one photon (q2 = q) is virtual and another one (q{) represents the external 
electromagnetic field and can be regarded as a real photon with the vanishingly small 
momentum q\ depends only on two invariant functions, longitudinal wl and transversal 
wt with respect to axial current index, 

T^ v x{qx,q) = ^ ~ W L (V) ^ 'qW T& 

+w T (V) (q 2 qiE pilv \ - <Ai<?2%«xA + <7 A <?i<?2%«™) • (13) 

Both structures are transversal with respect to vector current, q u T llu \ = 0. As for the 
axial current, the first structure is transversal with respect to q x while the second is lon- 
gitudinal and thus anomalous. 

In the local theory the one-loop result for the invariant functions wt and wl i^ 

i-ioop _ i-loop _ 2N c f 1 daa(l - a 



P = 2 Wrp~ = / — \ 2 ' 2 , (14) 

3 Jo a(l - a)q 2 + mj 

where the factor N c /3 is due to color number and electric charge. In the chiral limit, 
mj = 0, one gets the result for space-like momenta q (q 2 > 0) wl (q 2 ) = 2wt (q 2 ) = 
2/q 2 . 

The appearance of the longitudinal structure is the consequence of the axial Adler- 
Bell-Jackiw anomaljE^ESl For the nonsinglet axial current ^4( 3 ) there are no perturbative 



and nonperturbativeEH corrections to the axial anomaly and, as consequence, the 

invariant function remains intact when interaction with gluons is taken into account. 
Recently, it was shown that the relation 

w LT (q 2 ) = w L (q 2 ) - 2w T (q 2 ) = 0, (15) 
c Here and below the small effects of isospin violation is neglected, considering my = m 



which holds in the chiral limit at the one-loop level ©, gets no perturbative corrections 

from gluon exchanges in the iso-triplet case ^ . Nonperturbative nonrenormalization of 

RTI 

the nonsinglet longitudinal part follows from the 't Hooft consistency condition 1 ^, i.e. 
the exact quark-hadron duality realized as a correspondence between the infrared singu- 
larity of the quark triangle and the massless pion pole in terms of hadrons. OPE analysis 
indicates that at large q the leading nonperturbative power corrections to wt can only 
appear starting with terms ~ 1 /q 6 containing the matrix elements of the operators of di- 

mension six 1 ^. Thus, the transversal part of the triangle with a soft momentum in one of 
the vector currents has no perturbative corrections nevertheless it is modified nonpertur- 
batively. However, for the singlet axial current A^ due to the gluonic Ua (1) anomaly 
there is no massless state even in the chiral limit. Instead, the massive rj meson appears. 

So, one expects nonperturbative renormalization of the singlet anomalous amplitude Wj* 

at momenta below rf mass. Below we demonstrate how the anomalous structure w^ is 
saturated within the instanton liquid model. We also calculate the transversal invariant 
function wt at arbitrary space-like q and show that within the instanton model in the 
chiral limit at large q 2 all allowed by OPE power correct ions to wt cancel each other and 

only exponentially suppressed corrections remairJ^^SI The nonperturbative corrections 
to wt at large q have exponentially decreasing behavior related to the short distance 
properties of the instanton nonlocality in the QCD vacuum. 

Within the instanton liquid model the nondiagonal (VAV) correlator of vector current 
and nonsinglet axial-vector current in the external electromagnetic field is given by 



d A k 



Tr [T^k + q^Sik + q^ 



■T\ (k + qi ,k-q 2 )S(k- q 2 ) T v (k, k-q 2 )S (k) , 



(16) 



where the quark propagator S 1 (k) = k — M{k). The effective vector and axial-vector 
vertice s consist of the local and nonlocal parts and satisfy the Ward-Takahashi identities 
| 22 | 34 | -j^g s t ruc t ure f the vector vertices guarantees that the amplitude is transversal 
with respect to vector indices T au \(qi, q 2 )q\ = T^ v \(qi,q 2 )q 2 = and the Lorentz 
structure of the amplitude is given by (fT3l . 

The contribution of the diagram where all vertices are local to the invariant functions 
at space-like momentum transfer, q 2 = q 2 , are given by 
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(17) 



(18) 



where we also consider the combination of invariant functions wlt, (H31 . which show 
up nonperturbative dynamics. The notations used here and below are k + = k,k^ = 

k-q,k 2 ± = k + k_ - (fc+9 y~ g) . At large q 2 one gets wf c) (q 2 - oo) = showing 

saturation of the anomaly. The reason is that the leading large q 2 asymptotics of (flTt 
is given by the configuration where the large momentum is passing through all quark 



lines. Then the dynamical quark mass M(k) reduces to zero and the asymptotic limit 
of triangle diagram with dynamical quarks and local vertices coincides with the standard 
triangle amplitude with massless quarks and, thus, it is independent of the model. 

The contribution to the form factors when the nonlocal parts of the vector and axial- 
vector vertices are taken into account is given by 



w 



(nonloc) 




(19) 



Summing analytically the local ( flTt and nonlocal dT9l parts provides us with the result 
required by the axial anomaly ^ 



2 . 2N C 1 



3 q 



2 ' 



(20) 



Fig. [3] illustrates saturation of the anomaly in the non-singlet and singlet cases. Note, 
that at zero virtuality the saturation of anomaly follows from anomalous diagram of pion 
decay in two photons. This part is due to the triangle diagram involving nonlocal part 
of the axial vertex and local parts of the photon vertices. The result d20t is in agreement 
with the statement about absence of nonperturbative corrections to longitudinal invariant 
function following from the 't Hooft duality arguments. 

For wlt{<1 2 ) a number of cancellations takes place and the final result is quite simple 
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The behavior of u>lt(<? 2 ) is presented in Fig. |4] In the above expression the integrand 
is proportional to the product of nonlocal form factors / (k\) f (/c 2 ) depending on quark 
momenta passing through different quark lines. Then, it becomes evident that the large 
q 2 asymptotics of the integral is governed by the asymptotics of the nonlocal form factor 

/ (q 2 ) which is exponentially suppressed EES] Thus, within the instanton model the 
distinction between longitudinal and transversal parts is exponentially suppressed at large 
q 2 and all allowed by OPE power corrections are canceled each other. The instanton 
liquid model indicates that it may be possible that due to the anomaly the relation ( fT3t is 
violated at large q 2 only exponentially. 

The calculations of the singlet VAV correlator results in the following modification 
of the nonsinglet amplitudes 1 ^ 
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Figure 3: Normalized to l invariant function in the non- 
singlet case (solid line) and singlet case (dashed line). 



Figure 4: Normalized wlt invariant function versus Q 
predicted by the instanton model in the nonsinglet case 
(solid line) and singlet case (dashed line). 



Fig. |3] illustrates how the singlet longitudinal amplitude is renormalized at low 

momenta by the presence of the Ua (1) anomaly. The behavior of w^(q 2 ) is presented in 

Fig. 0] Precise form and even sign of w^(q 2 ) strongly depend on the ratio of singlet and 
nonsinglet couplings G'/G and has to be defined in the calculations with more realistic 
choice of model parameters. 



4 Conclusions 



We have analyzed the gluon field strength correlator (nonlocal gluon condensate), the 
vector Adler function and the nondiagonal vector - axial-vector correlator in external 
field (anomalous triangle amplitude) for Euclidean (spacelike) momenta within an effec- 
tive nonlocal chiral quark model motivated by the instanton model of QCD vacuum. The 
results on the gluon correlator are in good qualitative agreement with the lattice QCD 
simulations. The dominant contributions to the Adler function and the triangle ampli- 
tude come from the loop contributions of the light quark with dynamical momentum 
dependent mass. It is this contribution that provides the matching between low energy 
hadronized phase and the high energy QCD which is clearly seen in the behaviour of the 
Adler function (Fig. |2j and the wlt combination of the triangle amplitudes (Fig. @). The 
results obtained are close to estimates of the vector Adler function extracted directly from 
the ALEPH data on hadronic inclusive r decays and transformed by dispersion relations 
to the spacelike region. 
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